Bloch oscillations in an aperiodic one-dimensional potential 
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We study the dynamics of an electron subjected to a static uniform electric field within a one- 
dimensional tight-binding model with a slowly varying aperiodic potential. The unbiased model is 
known to support phases of localized and extended one-electron states separated by two mobility 
edges. We show that the electric field promotes sustained Bloch oscillations of an initial Gaussian 
wave packet whose amplitude reflects the band width of extended states. The frequency of these 
oscillations exhibit unique features, such as a sensitivity to the initial wave packet position and a 
multimode structure for weak fields, originating from the characteristics of the underlying aperiodic 
potential. 
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I. INTRODUCTION 

The nature of one-electron eigenstates has a significant 
influence on the electronic transport properties of solids. 
In pure periodic systems, the one-electron eigenstates are 
Bloch waves which are translational invariant and all de- 
localized in the thermodynamic limit. In the absence 
of scattering, the system behaves as a perfect conduc- 
tor whenever the Fermi energy falls into the conduction 
band. Disorder, originating from lattice imperfections, 
modifies the nature of the one-electron eigenstates. For a 
relatively weak disorder (of a magnitude smaller than the 
band width), the states at the band center may remain 
extended ?i but looses the phase coherence at large dis- 
tances. Those states, which lie near the band edges, turn 
out to be exponentially localized. Well defined energies, 
which separate localized and extended states, are known 
as mobility edges^ The metallic or insulating character of 
the system at zero temperature depends now on whether 
the Fermi level is located within the phase of delocalized 
or localized states, respectively. In the former case, the 
system shows a finite conductivity, while in the latter 
one it behaves as a perfect insulator. A large compared 
to the band width disorder localizes all the one-electron 
eigenstates, and the conductivity of the system vanishes 
at zero temperature. 

The above picture holds for three-dimensional (3D) 
systems (for an overview see Refs.|3 and0). In lower di- 
mensions the effect of disorder is much more dramatic. In 
particular, uncorrelated disorder of any magnitude causes 
exponential localization of all one-particle eigenstates in 
one dimension (1D)^& and weak localization in two di- 
mensions (2D) & At the end of eighties and beginning 
of nineties it was realized, however, that extended states 
may survive in ID systems when the disorder distribution 
is correlated ii2i2ii24iiii2iiiiii& Thus, a short-range cor- 
related disorder was found to cause the extended states 
at special resonance energies; they form a set of null 
measure in the density of states in the thermodynamic 



limit 
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Because of that mobility edges do not ex- 



ist. Oppositely, long-range correlations in the disorder 
distribution support a finite fraction of the delocalized 
statesi&ii and give rise to the existence of mobility edges. 
In Refs. Holl7lll8l it was argued that a nonrandom long- 
range inter-site coupling represents a one more driving 
force to delocalize one-particle states in iQiLiLiS and 
2Di£ geometries. The theoretical predictions about sup- 
pression of localization in ID geometry due to correla- 
tions in the disorder distribution were recently confirmed 
experimentally in semiconductor superlattices with in- 
tentional correlated disorder— as well as in single-mode 
wave guides with inserted correlated scatterers— 

Another class of ID models, that can exhibit an 
Anderson-like localization-delocalization transition, in- 
volves a nonrandom, deterministic potential which is in- 
commensurate with the underlying lattice. Several mod- 
els of this type have been extensively investigated in the 
literature, and the localized or extended nature of their 
eigenstates has been related to general characteristics of 
the incommensurate potentials— i2°*2ii22i22i An interest- 
ing class among them is represented by aperiodic slowly 
varying potentials. The latter support the mobility edges 
in a very close analogy with the standard 3D Anderson 
model— 

Recently, there appeared a renewed interest in the dy- 
namics of an electron in crystals subjected to a uniform 
static electric field. Under this condition, the electronic 
wave function displays so-called Bloch oscillations, 24,25,26 
the amplitude of which is proportional to the band width. 
The electronic Bloch oscillation was observed for the first 
time in semiconductor superlattices^ 7 - (for an overview 
see Ref . l2cf) . It is to be noticed that in bulk materials, 
this coherent regime of electron motion is hard to realize 
because the electron dephasing time due to scattering on 
defects and phonons is usually shorter than the period of 
Bloch oscillations. An exception represents semiconduc- 
tor superlattices, in which the opposite situation takes 
place. The unit cells in these materials are large enough 
to make the period of Bloch oscillations shorter than the 
electron dephasing time, so that several periods of os- 



dilations can be detected^ A similar phenomenon of 
sustained oscillations of the electromagnetic field, named 
photon Bloch oscillations, have also been reported in two- 
dimensional wave guide arrays and optical superlattices 
based on porous silicon^ 

Recently we investigated theoretically Bloch oscilla- 
tions in a ID disordered system with diagonal long-range 
correlated disorder and found that this type of correla- 
tions in disorder does not destroy the coherence of Bloch 
oscillations^ The amplitude of the latter was found to 
carry information about the energy difference between 
the two mobility edges. This result resembles the one 
that exists for an ideal Bloch band, where the amplitude 
of oscillation is proportional to the bandwidth. 

In this work, we further contribute to the general un- 
derstanding of the phenomenon of electronic Bloch oscil- 
lations in non-periodic low-dimensional systems exhibit- 
ing mobility edges. To this end, we will focus on the 
biased wave packet dynamics of a single electron moving 
in a lattice with an aperiodic slowly varying potential. 
As potential correlations are long-ranged, sustained os- 
cillations may be anticipated, similarly to our previous 
results^ It should be noticed, however, that the deter- 
ministic nature of the potential and its aperiodicity intro- 
duce new features to the Bloch oscillations. We solve nu- 
merically the time-dependent Schrodinger equation and 
compute the time dependence of the average electron po- 
sition. Fourier analysis is implicated to reveal regularities 
of the electron motion. 



II. MODEL SYSTEM AND RELEVANT 
MAGNITUDES 

We consider a tight-binding Hamiltonian on a regular 
ID open lattice of spacing a with an aperiodic slowly 
varying potential and a uniform static electric field ^&2i 

N 



TL = 2, \£n — eJ-an) \n)(n\ 

AT-1 

-jW|n)(n + l| + |n + l)<n|) , (i) 



where \n) is a Wannier state localized at site n with en- 
ergy e„, J- is the external uniform electric field and — e is 
the charge of the particle. The hopping amplitude is as- 
sumed to be uniform over the entire lattice with J > 0. 
In terms of the Wannier amplitudes ip n (t) = (n\^(t)), 
the time-dependent Schrodinger equation reads^& 



iipn = (e« - Fn)ip n - ip n 



Ipn 



(2) 



where we introduced the dimensionless magnitudes e n = 
s n /J and F — eJ-a/ J . Time is expressed in units of %/ J . 
Here, we will consider an on-site potential e n given 
by^ 



V cos ( -nan 1 ' j 



with V, a and v being variable parameters (an example 
of e n is depicted in Fig.QJ. For v = 1 this is just Harper's 
model, for which a rational a describes a crystalline solid, 
whereas an irrational a results in an incommensurate po- 
tential. Dynamical localization in the Harper's model 
under the action of electric fields was studied in Refl33. 
In the following, we restrict our study to an interesting 
range of parameters < v < 1 and < V < 2. It was 
demonstrated for this case that there exists a phase of 
extended states near the band center, which is separated 
by the two mobility edges ±E C = ±(2 — V) from two 
phases of localized states closer to the band edges^ 
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FIG. 1: (a) Aperiodic slowly varying potential calculated for 
V = 1, na — 1, and v — 0.5; b) The potential derivative 
showing the decrease of its local maxima on increasing the 
site index n. 

The pertinent quantities we will use to characterize 
the dynamics of the electron wave packet are its mean 
position (centroid) 



JV 

x{t) = Y(n 
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and the spread of the wave function (square root of the 
mean squared displacement) 
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V>0 , 
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where (n(t)) = YliLi n lV'n(' : )| 2 - As the initial packet is 
assumed spatially narrow, one has contributions to the 
wave packet dynamics, coming from a wide spectrum of 
eigenstates of the Hamiltonian 0J . 

In ideal lattices, a uniform field causes the electron 
wave packet to oscillate in space and time. The ampli- 
tude and the period of these oscillations are estimated 
semiclassically as L f = W/2F and tb = 2tt/F, respec- 
tively (see, e.g., Ref. l33|) . where W here is the width of the 
Bloch band in units of the hopping constant J. Subse- 
quently, the frequency of the harmonic motion is u = F, 
for the chosen units. The above picture was shown to 



remain valid when long-range correlated disorder is in- 
troduced with W being replaced by the band width W c 
of delocalized states^ In what follows, we provide nu- 
merical evidences of that the above semiclassical picture 
still holds for the present aperiodic system, although this 
requires a renormalization of the applied electric field due 
to the presence of the on-site aperiodic potential. 



III. NUMERICAL RESULTS 

We numerically solved the time-dependent Schrodinger 
equation @ by means of an implicit integration algo- 
rithm. 34 The initial wave packet was chosen to be a nar- 
row Gaussian of width a — 1, centered at an arbitrary 
lattice site no: 



^ n {0) = A(a) exp 



no) 



4ct 2 



(5) 



where A(o) is a normalization constant. To improve 
stability of the numerical algorithm, we set the electric 
potential to vanish at the initial site, replacing Fn by 
F(n — no) in Eq. (J2J. This produces only a shift in the 
origin of energy with no physical effect. 



A. Unbiased dynamics 

We start our analysis of the electron motion by study- 
ing the wave packet dynamics in the absence of the exter- 
nal field. The parameters of the aperiodic potential are 
set hereafter to V — 1, ira — 1, and v = 0.5. In Fig.|5fa), 
we depicted the time evolution of the centroid of a wave 
packet, which initially was located at site no = 48. It is 
worth to notice that the potential slope at this point is 
negative (see Fig.^J. The inset in Fig.|2{a) demonstrates 
that at the earlier stage of motion, the centroid oscillates 
harmonically, but loses the phase memory rapidly. The 
frequency of these few oscillations is about 0.13 in the 
dimensionless units. We show below that this number 
is directly related to the strength of the local field, pro- 
duced by the aperiodic potential in the vicinity of the 
initial position of the centroid, no = 48. Thus, the os- 
cillations found can be referred as to zero-field Bloch-like 
oscillations. To prove this statement, one should bear 
in mind that at short times the spread is small on the 
scale of the aperiodic potential [see Fig. Hfb)], so that 
the latter can be represented by two first terms of its 
Tailor expansion around the initial position of the cen- 
troid. Then, the potential derivative ds n /dn determines 
the local field strength which, in the dimensionless units, 
should be interpreted as the frequency of zero-field Bloch- 
like oscillations, similarly to the biased Bloch oscillations 
(see the preceding Section). Applying this reasonings to 
the site no = 48, we find that the potential derivative is 
about —0.14 at this site [see Fig. IHb)], whose absolute 
value is indeed fairly close to the frequency of oscilla- 
tions 0.13 extracted from the centroid dynamics. This is 
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FIG. 2: (a) Time-domain dynamics of the centroid of an unbi- 
ased wave packet (a — 1, no = 48 at t = 0). The inset shows 
the coherent oscillations at short times, (b) Time-domain dy- 
namics of the wave function spread for two different system 
sizes. 



an unambiguous confirmation of our qualitative picture. 
It is to be noticed, however, that these arguments fail 
in the vicinities of potential maxima and minima, where 
the potential slope is vanishingly small. 

As was already mentioned, the zero-field oscillations 
are not sustained for a long time. When the wave packet 
spread becomes comparable with the local scale of the 
aperiodic potential, the nature of the latter, with its wells 
and barriers, comes into play. This results in loosing the 
phase memory of the wave function. In Fig. [21b) we 
depicted the time-domain behavior of the wave packet 
spread. It displays a ballistic regime within some time 
interval presented by a linear in log-log scale dependence, 
followed further by a saturation at some value which, 
from one hand, depends on the system size, but from 
the other hand, is considerably smaller than that size. 
One can also see that the wave packet spread reveals a 
steep increase at some instant of time which is almost size 
independent. We relate this features to the tunneling of 
the wave packet from the potential well, where it was 
initially located, to the adjacent one. Figure |5] illustrates 
this tunneling process. 
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FIG. 3: Wave packet profile ^f n (t) (in arbitrary units) cal- 
culated for two moments of time, displaying the tunneling 
process from the initial potential well. All parameters are the 
same as in Fig. [5] 



B. Bias effects 
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FIG. 4: (a) Time-domain dynamics of the centroid of a biased 
wave packet (no = N/2 and a — 1 at t — 0) for two values of 
the applied electric field F — 1 and 1.25. (b) Fourier trans- 
form of the centroid. 



Once a uniform static electric field is turned on, the 
wave packet dynamics presents quite distinct features as 
compared to the unbiased behavior. In Fig.QJa) we plot- 
ted the centroids of wave packets calculated for distinct 
field strengths for the initial position at the chain center. 
It should be noticed that now Bloch oscillations remain 
sustained, i.e., no dephasing is taking place. Second, 
the oscillation amplitude is proportional to 1/ F as pre- 




8000 



FIG. 5: Time-domain dynamics of the wave function spread 
of a biased wave packet (no = N/2 and <j — 1 at t — 0). 
Notice its localized and breathing nature. 



dieted semiclassically. To provide further confirmation 
of the semiclassical picture, we calculated numerically 
the Fourier transform of the centroid, x(u>), as shown in 
Fig- Elk). Again, the estimated predominant frequency 
of the Bloch oscillations oj = F is corroborated. The lo- 
calized character of the wave packet is revealed in Fig. [S] 
where the spread of the wave function is shown; it satu- 
rates at a finite value independent of the chain size, dis- 
playing however periodic oscillations. The latter charac- 
terize the breathing nature of the oscillating wave packet. 

As already mentioned, the oscillations are not damped 
but remain amplitude-modulated. Solid lines in Fig-Ua) 
bound the spatial region within which the wave packet 
oscillates for a long time. The extent of this region L c is 
found to be L c ~ W c /F, where W c is independent of the 
applied field F. From the data in Fig. [3]we obtain W c ~ 
2E C . This value agrees remarkably well with the width 
of the band of extended states reported in Ref. 22. Thus, 
we arrive at one of the main conclusions of this work, 
namely there exist clear signatures of Bloch oscillations of 
a biased Gaussian wave packet between the two mobility 
edges. 

Some of the above characteristics of the Bloch oscilla- 
tions can depend on the initial location of wave packets. 
In particular, at the chain center the underlying aperiodic 
potential has an almost null derivative. Under this con- 
dition, the bias is solely due to the external field. In gen- 
eral, there is a contribution to the external field coming 
from the aperiodic potential; the total field at position tiq 
is given by F e g = de n /dn + F. The local contribution to 
the applied electric field shall be relevant whenever the 
potential gradient (in appropriate units) is of the order 
of the applied field. This is illustrated in Fig. HJa) with 
wave packets starting at locations with positive and neg- 
ative potential gradients. One can see that the oscillation 
amplitude becomes larger near regions with positive gra- 
dients and smaller for negative ones as compared with 
the results obtained in Fig. 0] for the same applied elec- 
tric field. Further, the characteristic frequency of these 
oscillations also exhibits shifts which are of the order of 
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FIG. 6: (a) Time-domain dynamics of the centroid of a biased 
wave packet (<r = 1 at t — 0) for F = 1 and two different initial 
positions, (b) Numerical Fourier transform of the centroid. 




FIG. 7: (a) Time-domain dynamics of the centroid of a biased 
wave packet (no = 48 and a = 1 at t = 0) and F — 0.5. 
(b) Numerical Fourier transform of the centroid. The multi- 
mode structure reveals the spread of the wave function on 
multiple potential wells. 



the local potential gradient [see Fig.^b)], as depicted in 

Fig. Eft). 

As a final remark, the Bloch oscillations shall develop 
a more complex structure for wave packets initially lo- 
cated near the beginning of the chain, especially at small 
electric fields for which the amplitude of the oscillations 
is large enough. Then the wave packet can probe regions 
with distinct scales of the potential gradient. In Fig.0we 
show a particular case of the dynamics of a wave packet 
starting at n — 48 under the applied field F = 0.5. 
Notice the rich pattern of the centroid oscillations. Its 



Fourier transform exhibits multi-mode structure with a 
series of narrow peaks which are shifted from u> — F by 
quantities that are of the order of the maximum poten- 
tial gradients at the region covered by the wave packet. 
The shifts to smaller frequencies can be much easily de- 
tected once the initial packet is concentrated at a region 
of negative potential gradient. However, the correspond- 
ing shifts to higher frequency can still be detected too as 
a series of small peaks. 



IV. SUMMARY AND CONCLUDING 
REMARKS 

We studied theoretically the single-electron wave 
packet dynamics within a tight-binding model with an 
aperiodic slowly varying site potential under the influ- 
ence of a static uniform electric field. The on-site poten- 
tial parameters were chosen such that, in the absence of 
the electric field, the model supports a phase of delocal- 
ized states at the center of the band, similarly to high 
dimensional disordered systems. The electric field pro- 
motes a bias which localizes the electron states. The re- 
sulting wave packet dynamics reveals Bloch oscillations. 
However, contrary to what occurs in disordered systems, 
where scattering on site potential fluctuations gradually 
degrades the oscillations, these remain sustained with no 
signature of depletion. 

The slowly varying aperiodic nature of the potential 
results in an additional local bias. Thus, the total bias 
has contributions of both the aperiodic potential and the 
one of the applied field. By defining an effective local 
field, we could show that the amplitude of the oscillations 
agrees well with the semiclassical prediction and can be 
used to estimate the width of the band of extended states. 
The typical frequency of these oscillations can also be 
understood on the semiclassical grounds. In particular, 
it is shifted depending on the local potential gradient. 
In the week field limit, these oscillations exhibit a multi- 
mode structure as the wave packet probes a larger region 
of the aperiodic potential. 

Our findings indicate that Bloch oscillations can indeed 
be observed in superlattices with slowly varying period- 
icity. The richness of the predicted dynamical behavior 
can lead to new electro-optical devices, aiming to explore 
the coherent motion of confined electrons. We hope that 
the present work will stimulate experimental activities 
along this direction. 
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